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Based on a realistic spin Hamiltonian for a frustrated quasi-two dimensional spin-f/2 antiferro- 
magnet Cs2CuCk, a three-dimensional spin ordering in the applied magnetic field B near the satu- 
ration value B c is studied within the magnon Bose-Einstein condensation (BEC) scenario. With the 
use of a hard-core boson formulation of the spin model, a strongly anysotropic magnon dispersion 
in CS2CUCI4 is calculated. In the dilute magnon limit near B c , the hard-core boson constraint is 
resulted in an effective magnon interaction which is treated in the Hartree-Fock approximation. 
The critical temperature T c is calculated as a function of a magnetic field B and compared with 
the phase boundary T c (B) experimentally determined in CS2CUCI4 [Phys. Rev. Lett. 95, 127202 
(2005)]. 



I. INTRODUCTION 

The insulator CS2CUCI4 is a quasi-two dimensional 
spin-1/2 antiferromagnet (AFM) composed of triangu- 
lar lattices (frc-plane) which are weakly coupled along 
crystallographic a-directiond In spite of a low spin, low 
dimensionality and geometrical frustration of spin inter- 
actions, a three-dimensional (3D) magnetic long-range 
ordering occurs in CS2CUCI4 at sufficiently low temper- 
atures T < Tjy — 0.6K. In the ordered state, magnetic 
moments lie in an easy 6c-plane and form an incommen- 
surate spiral structured An easy-plane spin anisotropy in 
CS2CUCI4 is produced by a Dzyaloshinskii-Moriya (DM) 
interaction with a DM vector T> || a. The anisotropy 
breaks SU(2) symmetry, however, U(l) symmetry cor- 
responding to spin rotations around T> vector is still 
present. The applied magnetic field B \\ a preserves U(l) 
symmetry and causes the spiral spin order to develop into 
a cone spin structure with increasing B: a longitudinal 
spin component grows with B, while a spiral transverse 
component decreases steadily and disappears completely 
at a saturation field B c . For higher fields, B > B c , the 
Zeeman interaction overcomes the antiferromagnetic spin 
coupling and the system enters a field-induced ferromag- 
netic stated In this state, ferromagnetic magnons become 
elementary spin exciatations with a gapped spectrum and 
a quadratic energy dispersion near magnon band minima. 
In the opposite direction, i.e. by starting from a high field 
B > B c and passing through B c , one finds a reversed be- 
havior: the field-induced magnon gap dissapears at B c 
and a spiral transverse magnetic ordering develops below 
B c . This ordered state is characterized by a gapless spin 
excitation spectrum with a linear low-energy dispersion. 

Such a phase transition observed^ 3 , in CS2CUCI4 in the 
applied magnetic field B || o in a vicinity of B c can be re- 
garded as a Bose-Einstein condensation (BEC) of a dilute 
gas of magnons. 4 '^ This suggestion was confirmed exper- 
imentally and supported theoretically in Ref. [jj. In this 
context, two particular features of CS2CUCI4 have to be 



mentioned. First, because of a rather weak dominant ex- 
change coupling J ~ 4.3K in this material, an exception- 
ally low saturation field B c ~ 8.5T can be easily achieved. 
Second, because of the extremely strong anisotropy of 
the ferromagnetic magnon spectrum in CS2CUCI4, 3D 
quadratic magnon dispersion is seen at very low energies 
E < E* ~ 50mK, and for E > E* the actual magnon 
spectrum is of 2D character. Therefore, an access to very 
low temperatures T ~ 50mK is required to be as close 
as possible to the asymptotic regime where universal 3D 
scaling laws are expected to hold. 

CS2CUCI4 is an easy-plane quantum AFM with a gap- 
less spin excitation spectrum at zero magnetic field. 1 Up 
to now, most of studies of a magnon (triplon) BEC have 
been reported 6 i 7 i 8 iV° on TlCuCl 3 and BaCuSi 2 6 that 
belong to a more rare class of gapped AFMs. In these ma- 
terials, a transverse spin long-range ordering occurs in the 
field range B cl < B < B c2 . In TlCuCl 3 and BaCuSi 2 6 , 
a saturation field B c i is rather high ^50T, which makes 
an observation of a magnon BEC at B c i a complicated 
matter, while lower fields near B c \ are accessible. Within 
the 3D magnon BEC scenario jiLidi 3 . near B c \ the or- 
dering temperature T c (B) is described by a power-law 
T c (B) ~ (B — Bc) 1 ^ with a universal critical exponent 
4>bec = 3/2. A few experimental findings show, how- 
ever, that the observed in TICUCI3 phase transition de- 
viates from a pure magnon BEC: first, an anisotropic 
spin coupling (of unknown nature) might break the axial 
U(l) symmetry of the system and produce a small but 
finite spin gapidi£*i& in the ordered state at B > B c \ 
and, second, the reported critical exponent <j) is some- 
what larger than predicted by theoryd Magnetic proper- 
ties of BaCuSi2C>6 are suggested to be well described by 
a quasi-two dimensional isotropic Heisenberg model and 
the reported^ recently convergence of the measured ex- 
ponent 4> to the predicted value 4>bec — 3/2 is in accord 
with the 3D magnon BEC scenario. 

In this paper a theoretical description of the magnon 
BEC transition in CS2CUCI4 near the saturation field is 
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presented with more details as compared to the previous 
work£ In section II, a model spin Hamiltonian applicable 
to CS2CUCI4 is discussed and rewritten in the hard-core 
boson representation. Two branches of the bare magnon 
dispersion and their densities of states are calculated. In 
section III, the hard-core boson constraint is treated in 
ladder approximation and the Bethe-Salpeter equation 
is solved to calculate the effective magnon interaction in 
the dilute gas limit. In the next section IV, the obtained 
effective interaction is considered in a mean-field approx- 
imation and the leading effects of magnon interaction are 
extracted and discussed. In section V, the critical tem- 
perature T c (B) as a function of applied magnetic field is 
calculated and compared both with the predicted univer- 
sal behavior of the phase boundaryii*i2iiii near B c and 
with experimental data. Short concluding remarks can 
be found in section VI. 



II. SPIN HAMILTONIAN AND MAGNON 
DISPERSION 

The spin Hamiltonian H of CS2CUCI4 involves the 
isotropic exchange H.q, the DM anisotropic term H.dm 
and the Zeeman energy TIb- By using the notations of 
Ref. IT?! the isotropic part Tto can be written as 

Ho = J^[>/5rSr +Ti+T2 + J'(Sr5r +Ti + 5rS*r +T2 ) 

+ J"S , rS , r +T3 ], (1) 



R 



where the nearest neighbor vectors t± and T2 are indi- 
cated in Fig. 1 and the out-of-plane vector T3 connects 
spins on vertical bonds between adjacent layers. A small 
relative shift of adjacent layers is neglected. Within each 




The crystallographic a-direction is chosen to be z axis 
and because V = (0, 0, V) the fee-plane becomes the easy 
plane for the ground state spin alignment in zero mag- 
netic field. In Eq. J5J a layer index n is introduced to 
show that DM vectors alternate between adjacent layers. 
In the present study we consider the applied magnetic 
field B || a and Zeeman energy reads 



-gfi B B 
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where g ~ 2.2 and [Ib is the Bohr magneton. The 
spin coupling parameters determined with high accuracy 
by neutron scattering measurements^ are: J = 4.34K, 
J'/J = 0.34, J" I J = 0.045 and V/J = 0.053. The total 
Hamiltonian Tt has U(l) symmetry arising from spin ro- 
tations around z axis. In CS2CUCI4 the saturation field 
B c ~ 8.5T and for B > B c at zero temperature the 
spins are fully polarized. To consider the magnetic phase 
transition in the vicinity of B c and at low T, we use 
the hard-core boson representation for spin- 1/2 opera- 
tors. Two types of bosons or and 6r/ are introduced for 
even and odd layers, respectively. This is given with a re- 
placement Sj[ — * br, — > a R and S R = 1/2 — o^or in 
even layers and the same transformation leads to 6r' for 
R/ belonging to odd layers. Such a replacement must be 
complemented by a constraint of no sites with boson oc- 
cupancies higher than unity. The constraint is satisfied 
by adding to TL an infinite on-site repulsion, IA — > 00, 
among the bosons 

H { u ] + H% ] = U ]T a+a+a R a R + U £ b+,b+,b w b w . 



r 



R' 



( 4 ) 

In the reciprocal space with an ortogonal basis we choose 
the Brillouin zone with wave-vectors restricted to ^ 
q x < 2-7T, ^ q v < 4-7T, and ^ q z < 2-7T, and define 
Fourier transforms of boson operators 



-iya R e' iqR , & q = -Ly&R, 



-zqR' 



R 



, (5) 



where Af = N/2 and N is the total number of the spin- 
1/2 lattice sites in the system. Fourier transforms of spin 
interactions now read 

J q = J cos q x + 2 J' cos (q x /2) cos (q y /2) , 



FIG. 1: Fragment of the triangular lattice (6c-plane) with 
nearest-neighbor exchange couplings, J and J' ~ J/3. The 
Dzyaloshinskii-Moriya coupling T> is on the same zig-zag 
bonds as the exchange J' 

layer, the dominant exchange J distinguishes a chain di- 
rection and the DM interaction is situated on interchain 
zig-zag bonds 

Hum - J2 ^ R x [5 R+T1 + Sr +T2 ]. (2) 

R 



J q ' = J"cos(g 2 /2), (6) 

£> q = 2£>sin(fe/2)cos (q y /2) . 

Two branches of 2D noninteracting magnons are de- 
scribed by the energy dispersions 

Eq" = ^q T Z>q + Afi, (7) 

where Ae = gy, B B - ( J + 2J' + J") . With these nota- 
tions, the bilinear part of the total Hamiltonian takes the 



3 



form 

q 

+ j£(o+ft q + 6+a q )] ) (8) 

and leads to two 3D magnon branches, A and B with 
bare dispersions 

4 B = J q T signPq^PqTjf + Ae. (9) 

The degenerate minima £q x = £q 2 are at Qi = 
(7T + <5i,0,0) for branch A and at Q 2 = (n ~ 8 2l 2ir, 0) 
for branch S. Without loosing precision we can use 
Si ~ 62 — 5 = 2 arcsin ( J'/2 J). The bilinear part of 
Ji now reads 

Hui =Y,[( E £- ^A+A^ + (E% - Mo )S+S q ] , (10) 
q 

where 



and 




In Eq. I|10fl. the bare chemical potential of magnons 
^0 = — [eqj + Ae] = — [ffq 2 + Ae] can be written as 
/io = .9Mb (-Be — B) , which defines a saturation field 
B c = W/ (gps)- Here the magnon bandwidth W is a 
function of spin coupling parameters J, J', J" and I?. In 
terms of J, W = 2.898J and therefore B c = 8.51T, as- 
suming g = 2.2. 

Due to the symmetry property of two magnon branches 
^Qi+q = ^Q2-q' tne i r densities of states (DOS) coinside, 
p-A \e) = p B (E), where 

q 

The upper end of the calculated DOS is at E u — W. For 
our purposes, the most important is the low-energy part 
of p A,B (E) depicted in Fig. 2. A critical point in the 
magnon spectrum is clearly seen at E — E* ~ 50mK: 3D 
character of the spectrum at E < E* changes into 2D 
one for E > E*. 
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FIG. 2: Low-energy part of magnon density of states (in units 
1/ J) for A and B branches. 

III. EFFECTIVE MAGNON INTERACTION 

When treating the magnon interactions, we take into 
account the dominant, hi — > 00, hard-core repulsion 
and neglect effects of a q-dependent part of interaction. 
Then, Him ~ H$ + , where 

o/(a) n,(b) _ 

n u u — 

( a qi a q 2 a q4 a q3 + ^qi^q 2 ^q4^q3)Aqi+q2,q3+q4i 

qi,q2, 
q3,q4 

(15) 

and A qi+q2j q 3+q4 implies the momentum conservation 

qi + q2 = q3 + q4- Below we refer to Tiffi and H$ 
as the hard-core magnon repulsion in a and b channels, 
respectively. 

In high magnetic field, B > B c , a fully spin polarized 
state is the vacuum state for magnons. Upon decreasing 
field, in the very vicinity of B c , i.e. B < B c and at low 
temperature the number of magnons is very small, n ~ 
(1 — B/B c ). In this case an effective magnon interaction 
can be found as a result of multiple magnon scattering 
by summing up ladder diagrams. 18 When accomplishing 
such a summation, we neglect interference between a and 
b channels. In this approximation, the problem reduces 
to solving the similar Bethe-Salpether equation for the 
renormalized scattering amplitudes r( a ) and r( h ) in both 
channels. Therefore, below we discuss with more detail 
the channel a. 

By substituting the transformation a q = a q .4 q — /3 q £> q 

into interaction Hamiltonian Ti^ in Eq. 11511 one finds 

that splits into sixteen scattering terms. In each 

process, the bare interaction acquires now an extra fac- 
tor given by a product of four a q and /3 q coefficients. For 
instance, a scattering process (Aq 3 , S q4 ) — > (-4 qi , £> q2 ) 
is described by the amplitude 2a qi /? q2 o; q 3/? q4 W. Taking 
into account multiple magnon scattering in the ladder ap- 
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proximation leads to a replacement of 21A by an effective 
two-particle interaction r^ a ^ (K, cj), where K = qi+q2 = 
Q3 + Q4 and oj — lui +L02 — LU3 + un are the total momen- 
tum and energy of a scattering magnon pair; the extra 
factor remains unchanged. The same holds for the other 
fifteen scattering processes. In the magnetic field close 
to B c and at low T only the particles near the magnon 
band minima at q = Qi,2 are present, which allows us to 
send w -> and consider r( a )(K, 0) = T 



(a) 
K • 



94 



'12 



r(°) I ( b ) 
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minima wave- vectors q = Qi, Q2 are occupied and con- 
tribute to the phase transition. In this section, we con- 
sider first how the magnon spectrum near the band min- 
ima is renormalized due to the effective magnon inter- 
action. The study is based on the Hartree-Fock (HF) 
approximation, since near the quantum critical point, 
B = B c and T = 0, the HF theory is suggested to de- 
scribe correctly the phase boundary>i2ii^ 

We assume that the temperature approaches T c from 
above. Then in the disordered phase, the mean-field in- 
teraction Hamiltonian in the channel a can be written 



as 



a/a) r (a) 

rL int,MF — J^J- 1 qi + 

qi.q2 



12 ( a q>qi) a <! 



(17) 



FIG. 3: Bethe-Salpeter equation for the scattering vertices 

p(a),(6) 

The effective interaction F^ is obtained in an ex- 
plicit form by solving the Bethe-Salpeter equation pre- 
sented in Fig. 3 diagrammatically. In this presentation, 
Qi = (qi: w i) and k — (k, a/) are four-dimensional vec- 
tors. Any external incoming or outgoing line is associated 
with either A or B magnon; two internal lines labelled 
by 94 + k and q% — k correspond to Greens's function 
(q±k) = a 2 q±k G^ {q±k)+ P 2 q±k Q^ (q±k). Here 
(q ± k) and (q ± k) are the Green's functions of 
decoupled A and B magnons. Since the bare interaction 
U is q- and ^-independent, the solution of the Bethe- 
Salpeter equation in the zero frequency limit is 



qi+q2^ 



<+k<-k 



a qi+k^q 2 -k 



B 2 a 2 

, ^qi+k a q 2 -k 



£q 2 -k 



^2-k 



^qi+k/^q 2 ~k 



^qi+k 



q2-k 



^qi+k 



q2-k 



.(16) 



The similar derivation applies also to an effective magnon 
interaction r q ^ +q2 in the channel b, which can be ob- 
tained merely from r q "l|_ q2 by interchanging A <-> B. 

By collecting back all the terms in both a and b chan- 
nels, one arrives at an effective interaction Hamiltoni- 
ans of magnons. It becomes clear now that the effective 
Hamiltonian takes the original form (|15|l where the bare 

parameter 2U has to be replaced by r q ^ + and r qi + 
in channel a and &, respectively. 



IV. MEAN-FIELD THEORY 

The main goal of the present study is to calculate a 
shape of the phase boundary, i.e. a field dependence 
of the critical temperature T c (B) for B < B c . When 
approaching B c from below, T c — > 0. In the low-T re- 
gion near £> c , only low-energy magnon states at the band 



which is approximated further as follows 

- aqi/^ltU «£ qi + B+Au) • (18) 



Ql ^ r S+m,0L r Q2+q2/ 3 Q2> SmCe l0W " 



(a) 



In this expression, we made replacements r q "^. 



p(a) 

energy magnon densities (?i q ) = (*4 q *4 q ) and (?i q ) 
(S q B q ) are located in the q-space near Qi and Q 2 for 
A and B branch, respectively. Because of a symmetry of 
the magnon interaction, the band minima degeneracy is 
preserved and one has jj £ q (n*) = ^ E q (™q) = n - 
A possible shift due to magnon interaction of the band 
minima location in the q-space at Qi,2 is assumed to bc 
small and neglected below. In the zero-order approxima- 
tion , the quantitie s (A ^ B cll + $ qi „4 qi ) vanish, and we 
omit them in Eq. (|18fl as well. A validity of the latter 
approximation will be verified later on. 

Similar arguments as above allow us to approximate 
the operator part in Eq. I|17|) as follows 



-a q /3 q (A+B^ - 



jB q S q 



(19) 



and to treat the effective magnon interaction in the chan- 
nel b in the same manner as in the channel a. Finally, 
by noting that a 2 Qi = a 2 Q2 = a 2 , 8^ = B 2 Q2 = B 2 and 
a 2 3> (3 2 , we retain only the leading terms in the mean- 
field interaction Hamiltonian, which leads to the follow- 
ing result 



H 



(a) 



int,MF ' rL int,MF ~ 

2r ? ^(^ q + 6+6 q 



-2n^r q (^+S q + 6+^l q ), (20) 
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where 

r = = a-T^ 

(21) 

r' = a 8 \v {a) - r (b) 1 

1 q ~ "qPql 1 Q 1+q i Q 2 +qJ 

According to Eq. I|20(l . the bare chemical potential of 
magnons is renormalized 

fi -> fi ef f = Mo - 2r?i. (22) 

With the use of the definition 112|) for a q and l|16l) for 
we obtain numerically the estimate T ~ 0.85 J. The 
second term in Eq. (20) generates magnon hybridization. 
This term shifts the bottom of the magnon bond slightly 
down and leads to a weak mass enhancment of low-energy 
magnons. Both effects are proportional to n 2 and we 
omit them since n -C 1 near T c for (B c — B) -C B c . 
In this region, the low-energy magnon A and B branches 
remain decoupled, which justifies an approximation made 
before. 

V. CRITICAL TEMPERATURE T c (B) NEAR B c 

In the previous sections we obtained in the dilute 
magnon limit n <C 1, and within the HF approximation 
that the leading, linear in n, effect of magnon interac- 
tion is a renormalization of the chemical potential, Eq. 
1)22(1. In this section we calculate the critical tempera- 
ture T c (B) of a magnon BEC transition as a function 
of B assuming that for given B < B c the temperature 
approaches the phase boundary from the normal phase. 




6.75 7.25 7.75 8.25 

B, T 

FIG. 4: Calculated phase boundary T c (B) near B c (solid 
line) . Experimental data (squares) are from specific heat mea- 
surements reported in Ref. y 

As follows from Ea. H22|l . a comparatively large magnon 
density n at high temperature makes the chemical poten- 
tial to be negative, /U e // < 0, and all the magnons are 
noncondensed. Under decreasing T, (J- e ff tends to zero 
and a BEC transition is expected to occur when fJ- e ff — 



at T c (B); for T ^ T c , a condensate density no starts to 
develop continuously at the magnon band minima. From 
Eq. (|22|l . the condition for the phase transition is written 
asi 

9[ ib (B c — B) = 2Tn (T c ) . (23) 

Here n (T) = A^ 1 £ q /s (£ q ) with f B (£ q ) being 
the Bose distribution function taken at ^ e // = and 
= E£ or = E®. This means that for T < T c 
the magnon condensate develops simultaneously at 
ct = Qi,2- Since at fi e ff = the distribution function 
Jb {E) diverges as T/E for E — > 0, a magnon spectrum 
near the band minima mainly contribute to n (T c ). If T c 
is sufficiently low, the very low-energy magnons having 
a 3D quadratic dispersion dominate and drive the BEC 
transition. In this asymptotic regime the power-law is 
expected 11 ! 12 ! 13 T c (B) - (B c - B) 1 ^ with the universal 
exponent 4>bec — 3/2. To check this suggestion and 
to compare the calculated phase boundary with that 
determined experimentally in CS2CUCI4, we solved the 
equation i|23|l with the use of the magnon density of 
states, Fig. 2, realistic for CS2CUCI4. The result of 
calculations, with B c — 8.51T and T — 0.85J, is shown 
in Fig. 4 and compared with experimental data— for 
B < B c . In the very vicinity to B Cl the theoretical curve 
is well reproduced by the power-law with the calculated 
(j>th — 1-5, while a fit of experimental data yields 3 - 
4>exp = 1.52(10); both values of the critical exponent are 
close to the universal 4>bec — 3/2. 



VI. CONCLUSIONS 

Based on a realistic quantum spin model in the hard- 
core boson representation, we developed a mean-field 
theory of a magnon BEC transition in CS2CUCI4 near 
the saturation field B c . The calculated critical temper- 
ature T c (B) describes surprisingly well the experimen- 
tally measured phase boundary in a narrow field region 
B c — AB < B < B c , where AB/B C ~ 0.06. In this re- 
gion, where ^ T c < 300mK, both the measured and 
calculated critical exponents (f> are very close to the uni- 
versal value 4>bec = 3/2 characteristic for 3D quadratic 
magnon dispersion. It means that in CS2CUCI4 the 
magnon BEC near B c is driven by 3D magnon states with 
energies below E* w 50mK. Actually, by sending for- 
mally E* — > (see Fig. 2) in our calculations we obtained 
T c = at any B. At lower fields B < B c - AB ~ 8T, a 
mean-field description of the experimentally determined 
phase boudary T c (B) fails, as was expected. Both ther- 
mal fluctuations and a temperature renormalization of 
the magnon spectrum have to be taking into account to 
improve a theoretical description of the phase boundary 
at lower magnetic fields. 
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